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PREFACE 
This  final  report  presents  an  executive  summary  of  the  three  interim 
reports  prepared  for  the  project  "Use  of  Synthetic  Demand  Modelling  Techniques 
in  Transportation  Planning  for  Small  Urban  Areas  in  Indiana".   The  titles  of 
these  three  reports  are: 

I.  "Framework  for  the  Transferability  of  Trip  Generation  Parameters 
of  Small  Urban  Areas" 

II.  "Evaluation  and  Development  of  Synthetic  Trip  Distribution  Modelling 
Techniques  as  Applied  to  Small  Urban  Areas  in  Indiana" 

III.  "Effects  of  Simplifying  Traffic-Zone  and  Street-Network  Systems  on 
the  Accuracy  of  Traffic  Assignments  in  Small  Urban  Areas  in  Indiana" 

This  summary  will  serve  as  a  set  of  guidelines  for  the  application  of 
synthetic  demand  models  in  small  urban  areas.   Recommendations  have  also 
been  made  as  to  the  use  of  these  guidelines. 

The  graduate  research  assistants  who  worked  on  this  project  were:   Hani 
S.  Mahmassani,  James  R.  Mekemson,  and  Edward  W.  Hanscom.   Another  graduate 
research  assistant,  Karen  L.  Pickett,  assisted  in  the  preparation  of  the 
final  report.   The  following  planning  division  staff  of  the  Indiana 
State  Highway  Commission  provided  much  help  and  assistance:   Mike  Faris,  John 
Mick  and  James  Barr.   The  typing  of  the  reports  was  most  ably  performed  by 
Marian  J.  Sipes  and  Melody  Gray.   The  constant  guidance  of  Professor  Harold 
L.  Michael  is  gratefully  acknowledged. 


INTRODUCTION 

Initiating  a  continuing,  comprehensive,  and  coordinated  transportation 
process  is  a  considerable  task  that  faces  small  metropolitan  areas  as  their 
population  reaches  50,000.   Full-scale  comprehensive  transportation  studies 
used  in  large  areas  have  technical  and  financial  requirements  that  often  ex- 
ceed the  capabilities  of  these  smaller  areas.   In  most  cases,  such  sophistica- 
tion is  not  needed  because  the  problems  faced  by  smaller  areas  are  different 
in  nature  and  in  complexity  from  those  in  larger  areas.   It  is  therefore 
imperative  to  develop  and  implement  simplified  alternative  planning  approaches 
which  would  greatly  reduce  the  time  and  cost  currently  required  by  the  "classical" 
transportation  planning  process.   One  approach  to  simplifying  the  transportation 
planning  process  is  to  eliminate  the  need  for  a  full  scale  origin-destination 
(0-D)  survey.   This  can  be  achieved  by  use  of  a  synthetic  travel  demand 
modelling  approach,  whereby  information  traditionally  obtained  from  the  0-D 
survey  is  fabricated  or  synthesized  through  the  use  of  parameters  "borrowed" 
from  "similar"  areas. 

The  first  phase  of  this  study  was  aimed  at  reducing  the  cost  and  complexity 
of  the  conventional  planning  procedures  by  developing  a  framework  for  trans- 
ferring trip  frequency  parameters,  thus  eliminating  the  need  for  conducting 
full-scale  travel  surveys.   To  this  end,  the  objectives  were: 

1.  Compilation  and  interpretation  of  various  travel  characteristics  for 
the  several  small  metropolitan  areas  in  Indiana  in  which  0-D  studies 
have  been  conducted. 

2.  Comparison  of  these  travel  characteristics  and  identification  of  the 
other  urban  variables  that  explain  these  characteristics. 

3.  Detailed  examination  of  trip  generation  parameters  and  models  for 
those  areas  at  three  levels  of  aggregation:   areawide,  zonal  and 
household;  assessment  of  the  stability  and  transferability  of  these 
parameters  for  synthetic  modeling  applications. 


4.   Development  of  a  generalized  framework  for  the  transferability  of 
trip  frequency  (generation)  parameters  and  models. 

The  second  phase  of  the  study  concentrated  solely  on  the  trip  distribution 
part  of  the  modelling  process.   The  most  commonly  used  trip  distribution  model 
is  the  gravity  model  in  which  calibration  of  the  model  primarily  involves  the 
establishment  of  the  appropriate  travel  time  impedance  function.   This  phase 
of  the  research  evaluated  the  performance  of  some  of  the  various  suggested 
synthetic  gravity  model  calibration  techniques  and  developed  an  improved  pro- 
cedure for  synthetic  self  calibration  of  the  gravity  model.   The  same  small 
urban  areas  of  Indiana,  as  were  examined  in  the  trip  generation  phase,  were  used 
here.   In  addition,  the  complete  transportation  modelling  data  sets  from  three 
of  these  urban  areas  were  used  for  detailed  analysis  and  testing.   The  five 
major  tasks  were  to: 

1.  Summarize  the  travel  time  statistics  of  the  study  areas  and  relate 
these  statistics  to  other  urban  area  characteristics.   Determine 

if  these  statistics  will  be  of  any  use  in  other  urban  areas  in  which 
the  transportation  planning  process  is  to  be  implemented  or  updated. 

2.  Investigate  the  sensitivity  of  the  gravity  model  with  respect  to 
friction  factor  errors  introduced  by  the  use  of  borrowed  or  standardized 
friction  factors. 

3.  Investigate  the  possibility  of  using  areawide  travel  time  statistics 
and  other  variables  in  a  self-calibrating  gravity  model. 

4.  Investigate  the  possibility  of  using  origin  zone  specific  travel  time 
statistics  in  a  synthetically  self-calibrating  gravity  model. 

5.  Evaluate  and  make  recommendations  as  to  what  method  of  synthetic 
trip  distribution  modelling  should  be  used  in  the  transportation 
planning  process  for  small  urban  areas. 


In  the  third  and  final  phase  of  this  endeavor,  a  study  was  made  to  deter- 
mine if  the  complexity  of  urban  travel  forecasting  could  be  reduced  by 
simplifying  traffic-zone  and  street-network  systems  without  sacrificing  accuracy 
in  traffic  assignments.   This  phase  of  the  study  considered  simplification  of 
several  forms  using  transportation  network  and  travel  data  from  two  small 
Indiana  urban  areas.  Alternative  traffic-zone  and  street  network  systems  were 
developed  and  tested  for  each  city,  then  test  results  were  compared  by  overall 
measures  of  traffic  distribution  and  assignment  accuracy,  and  by  statistical 
analyses  of  network  links  stratified  by  volume  group. 

In  the  following  sections  are  briefly  summarized  the  major  findings  of  the 
study  including  a  set  of  guidelines  for  the  application  of  synthetic  demand 
models  for  estimating  travel  demand  in  small  urban  areas  in  Indiana.   The  guide- 
lines indicate  the  recommendations  for  new  transportation  studies  or  updating 
existing  studies. 

SYNTHETIC  TRIP  GENERATION 

The  need  for  simplification  of  the  transportation  planning  process  in 
small  urban  areas  is  apparent.   Various  efforts  have  been  undertaken  in  this 
direction.   Recent  efforts  have  examined  the  use  of  cross  classification 
disaggregate  models  for  generation  analysis.   This  technique  has  been  strongly 
suggested  for  synthetic  modelling  applications,  especially  in  small  urban  areas. 
It  consists  of  developing  household  trip  rates  for  each  household  socio-economic 
category.   Since  the  household  is  taken  to  be  the  basic  unit  of  analysis,  these 
models  are  behavioral  in  nature  and  should  therefore  be  transferable  from  one 
area  to  another. 


The  stability  and  transferability  of  household  trip  rates  between  small 
and  medium  areas  is  one  of  the  matters  examined.   In  contrast  to  previous  work 
done  in  this  area,  this  study  developed  a  framework  for  transferring  travel 
frequency  parameters  at  three  levels  of  aggregation:   areawide,  zonal,  and 
household.   This  theoretical  framework  is  supported  by  empirical  evidence 
derived  from  the  comparison  of  parameters  obtained  from  various  study  areas  for 
each  of  the  relevant  levels  of  aggregation. 

General  Characteristics  Of  The  Study  Areas 
This  study  was  aimed  at  compiling  information  from  those  small  and  medium 
sized  metropolitan  areas  in  Indiana  in  which  a  detailed  travel  survey  has  been 
conducted  and  a  transportation  study  performed.   Seven  metropolitan  areas  with 
population  in  the  range  of  50,000  to  250,000  were  thus  included  in  this  study. 
These  are:   Anderson,  Evansville,  Fort  Wayne,  Lafayette,  Muncie,  South  Bend, 
and  Terre  Haute.   These  areas  are  Standard  Metropolitan  Statistical  Areas 
(SMSA's)  as  defined  by  the  U.S.  Bureau  of  Census,  and  the  general  characteristics 
of  these  urban  areas  are  shown  in  Table  1 . 

Distribution  Of  Trips  By  Purpose 
Trips  are  usually  stratified  into  three  types  in  urban  transportation 
planning  studies.   These  are: 

-  internal  trips  which  are  those  with  both  ends  within  the  study  area 
cordon  line, 

-  internal-external  trips,  with  one  end  inside  the  study  area  and  the  other 
outside  the  cordon,  and 

-  external-external  or  through  trips,  with  both  ends  outside  the  study 
area. 


9 
•a: 


en  oo 

•—       cm  i—       oo       o       oo 

r»        r^  ^r        •—        o        ^- 
cr>          « 

p-  CM  p—  r-  p-  CM 


Q 
Z 
UJ 
CD 


O 
00 


oo 

■— 

r-» 

r— 

s 

O 

in 

r-~ 

10 

o 

CO 

cm 

,-- 

CTi 

n 

n 

• 

■ 

• 

f-» 

o% 

10 

'"" 

■"" 

CO 

CM 


03 
O 

•H 
•U 
W 

•H 
p-i 

cu 

4J 

O 
« 
U 
CO 
J= 
c_> 


o 

z 


to 

r—  10 

r».       o 

en         « 

i—        o 

o 


IT) 

o 


o 

CO 


CM 

in 


CO 
CO 


CO 


>- 
< 


IT) 

CO 

o 

CM 

in 

P~ 

*3- 

o 

r» 

p— 

r>» 

CM 

in 

*r 

CTl 

M 

n 

• 

• 

• 

"~ 

o 

CTl 
CM 

"~ 

•"" 

CO 

ITS 
0J 


o 


o 
o 


CO 

c 


0) 

u 
< 

3 


>- 
< 


tx. 
o 


z 
< 


r- 

O 

p— 

te 

00 

10 

CO 

r^ 

Ov 

■1 

« 

pa- 

CO 

^, 

o 

10 

CM 

in        in        in 

CO  p—  r— 


p—  i—  CO 


*r 

O 

o 

r— 

O 

10 

o 

o 

r^ 

in 

in 

C\J 

CO 

CTl 

en 

■"" 

in 

o 

'■" 

l— 

CM 

d) 

z 

■H 

o 

XI 

00 

r* 

CO 

CO 

cfl 

eg 

r>. 

CO 

O 

o 

o 

CO 

H 

UJ 

en 

CO 

CO 

CO 

CO 

o 

o 

#— 

n 

* 

• 

• 

• 

z 
< 

O 

cn 

cn 

CM 

r— 

r~ 

CO 

OJ 

s_ 

<: 

>> 
-o 

3 
•!-> 
00 

s_ 

o 

o       u_ 

r»«. 

CTl  </> 

p—  — 

I  s- 

o 

<  Q. 

oo        oj 

S  Ce: 

oo 

0)        -o 


•r— 

■o 

c 

* 

p— 

=3 

* 

* 

o 

■o 

TS 

JZ 

o> 

p— 

p— 

0) 

c 

o 

o 

l/l 

•p— 

JE 

J= 

3 

f— 

OJ 

OJ 

o 

p— 

10 

01 

3= 

01 

3 

3 

s- 

$ 

o 

O 

i- 

(O 

c 

o 

3: 

3: 

OJ 

oj 

o 

o. 

^ 

■y— 

-o 

s_ 

(. 

-t-> 

OJ 

OJ 

OJ 

1/) 

>» 

IO 

•r— 

a. 

o. 

c 

oj 

p— 

O. 

o 

> 

3 

3 

o 

o 

in 

!_ 

o. 

u 

■4-' 

4-> 

s_ 

3 

o 

o 

3 

3 

OJ 

00 

a. 

o 

< 

< 

a. 

c 
o 


13  ■>-> 

—  S- 

o 

oo  — 

3  i/> 

i/i  C 

C  ITS 

0J  S- 

O  I— 

e  e 

o  o 


On  the  average,  82.2  percent  of  the  vehicle  trips  were  found  to  be  internal, 
15.2  percent  external-internal,  and  2.6  percent  through  trips.   These  percentages 
are  the  (unweighted)  averages  of  the  respective  percentages  for  each  of  the 
seven  urban  areas.   The  home  interview  survey,  the  most  costly  and  time  consuming 
element  of  the  0-D  survey,  is  the  basis  for  determining  internal  trips  for  the 
base  year.   Inasmuch  as  the  aim  of  this  research  effort  was  to  examine  synthetic 
techniques  which  would  minimize  or  eliminate  the  need  for  the  home  interview 
survey,  it  essentially  focused  on  internal  trips. 

Number  and  Definition  of  Trip  Purposes.   For  smaller  urban  areas,  stratifi- 
cation of  internal  trips  into  three  purposes  is  usually  adequate.   These  are: 
Home-Based  Work  (HBW) ,  Home-Based  Other  (HBO) ,  and  Non-Home  Based  (NHB) .   In 
Table  2  are  given  the  total  internal  vehicle  trips  by  purpose  along  with  their 
percentages.   It  was  observed  that  the  distribution  by  purpose  of  internal 
vehicle  trips  was  not  the  same  for  all  urban  areas  in  Indiana.   However,  for  an 
initial  assumption  (to  be  adjusted  later  in  the  process)  in  a  synthetic 
modelling  effort  the  following  unweighted  vehicle  trip  averages  can  be  made: 
HBW  -  16.5%;  HBO  -  46.9%;  NHB  -  36.0%. 

Variation  Within  Urban  Areas.   A  household  level  of  analysis  was  used  to 
investigate  some  of  the  factors  behind  the  variation  of  trip  purpose  distribution 
between  urban  areas.   A  Chi-square  test  was  used  to  test  the  hypothesis  that 
the  distribution  by  purpose  of  trips  made  by  households  was  independent  of  the 
socio-economic  characteristics  of  households  (as  reflected  in  auto  ownership 
status),  and  the  results  led  to  the  rejection  of  the  independence  hypothesis 
with  99.5%  confidence.   Therefore,  the  distribution  of  vehicle  trips  by  purpose 
turned  out  to  be  significantly  different  between  auto  ownership  groups. 

Inasmuch  as  the  trip  distribution  by  purpose  is  significantly  different  for 
each  socio-economic  group,  it  would  appear  reasonable  to  assume  that  the  overall 
trip  distribution  by  purpose  in  two  different  urban  areas  would  be  the  same, 


Table  2.   Distribution  of  Total  Vehicle  Trips  by  Purpose 


City 

HBW 

HBO 

NHB* 

TOTAL 

ANDERSON 

47,105 
(16.2) 

145,958 
(50.3) 

97,139 
(33.5) 

290,202 
(100) 

EVANSVILLE 

74,586 
(16.5) 

188,285 
(41.6) 

189,923 
(41.9) 

452,794 

(100) 

FORT  WAYNE 

90,203 
(23.3) 

146,990 
(38.0) 

150,029 
(38.7) 

387,222 
(100) 

LAFAYETTE 

44  ,337 
(15.0) 

154,067 
(52.1) 

97,306 
(32.9) 

295,710 
(100) 

MUNCIE 

38,591 
(13.4) 

159,017 
(55.2) 

90  ,641 
(31.4) 

288,249 
(100) 

SOUTH  BEND 

79,672 
(15.1) 

208,452 
(39.3) 

241,896 
(45.6) 

530,020 
(100) 

TERRE  HAUTE 

36,745 
(15.7) 

121,836 
(52.1) 

75,130 
(32.2) 

233,711 
(100) 

including  Truck  trips 

(Numbers  in  parentheses  are  percentages) 


Unweighted  Average  Percentages:  HBW:  16.5% 

HBO:  46.9% 
NHB:  36.6% 
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if  the  socio-economic  mix  in  these  areas  were  similar.   This  implies  that 
"borrowing"  areawide  percentages  without  first  comparing  the  socio-economic 
mix  of  the  two  areas  could  lead  to  erroneous  results. 

Areawide  Trip  Frequency 
"Areawide",  as  it  is  used  in  this  context,  describes  the  geographic  unit 
of  analysis  which  is  the  totality  of  the  urban  area  for  which  a  transportation 
study  is  being  conducted.   Areawide  frequency  measures  are  to  be  deduced  based 
on  aggregate  areawide  characteristics  of  the  urban  area.   The  following  rates 
are  the  most  commonly  used  descriptive  measures  of  trip  frequency: 

-  average  number  of  trips  per  household 

-  average  number  of  trips  per  person 

-  average  number  of  trips  per  automobile. 

The  vehicle  trips  per  household  ranged  from  5.64  to  9.94  per  day,  while  vehicle 
trips  per  person  varied  from  1.90  to  3.22.   This  variability  is  to  be  expected 
since  areawide  travel  demand  reflects  various  urban  area  characteristics  which 
are  most  likely  different  between  the  study  areas. 

Areawide  household  trip  frequency  parameters  depend  on  the  socio-economic 
mix  of  the  households  in  the  urban  area  under  study.   Thus,  in  order  to  transfer 
aggregate  areawide  parameters  from  one  area  to  another,  care  must  be  taken  to 
thoroughly  compare  the  socio-economic  characteristics  of  the  households  of  those 
two  areas.   For  gross  areawide  estimations,  transferring  aggregate  trip  frequency 
parameters  between  urban  areas  with  roughly  similar  socio-economic  characteristics 
is  feasible.   More  detailed  analysis  and  forecasting  at  a  finer  level  of 
aggregation  would  most  definitely  not  rely  on  borrowed  aggregate  areawide  travel 
frequency  parameters,  but  should  try  to  borrow  parameters  developed  at  a  finer 
level  of  aggregation,  such  as  the  traffic  zone  or  the  individual  household. 


Trip  Generation  Analysis  At  The  Zonal  Level 
Multiple  regression  techniques  are  generally  used  to  relate  average  zonal 
trips  to  the  zonal  socio-economic  and  land  use  characteristics.  All  seven  study 
areas  developed  zonal  trip  generation  equations.   The  trip  purposes  for  which 
these  equations  were  developed  varied  between  the  study  areas;  however,  these 
equations  were  combined  according  to  the  three  purposes  defined  for  this  study 
(HBW,  HBO,  and  NHB)  in  order  to  have  a  basis  for  comparison. 

The  transferability  of  the  internal  vehicle  trip  production  equations  were 
investigated  at  two  levels  of  detail: 

1.  For  each  urban  area,  the  mean  zonal  value  of  trips  produced  for  each 
purpose  were  predicted  using  the  equations  calibrated  in  other  areas. 
Areawide  mean  zonal  values  of  the  independent  variables  were  used  for 
these  predictions,  which  were  in  turn,  compared  to  the  corresponding 
observed  values  (from  the  0-D  survey) .   The  percent  difference  of  the 
mean  (difference  between  mean  of  estimated  zonal  trips  and  mean  of 
observed  zonal  trips  as  a  percent  of  the  observed  mean)  was  then 
calculated  as  a  measure  of  the  overall  predictive  ability  of  each 
equation  in  each  urban  area. 

2.  A  more  detailed  evaluation  was  made  for  selected  urban  areas  (based 

on  data  availability)  by  predicting  zonal  trips  for  each  purpose  using 
zonal  values  of  the  independent  variables.   These  predictions  were 
compared  on  a  zone-by-zone  basis  to  the  corresponding  0-D  observations. 
The  percent  standard  error  of  the  estimate,  expressed  as  a  percentage 
of  the  mean  value,  was  used  to  measure  the  overall  predictive  ability 
of  an  equation. 
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The  investigation  of  the  transferability  of  the  trip  attraction  regression 
equations  was  limited  to  comparisons  of  the  mean  zonal  trips  (observed  vs. 
estimated)  for  three  of  the  study  areas  (Lafayette,  Fort  Wayne,  Muncie) .   Over- 
all, trip  attraction  equations  can  be  transferred  for  the  purpose  of  developing 
rough  estimates  that  would  be  redefined  later  on  in  the  process. 

The  results  obtained  in  this  study  showed  that  zonal  aggregate  model 
parameters  tend  to  be  dependent  on  the  particular  zonal  scheme  for  which  the 
model  was  developed.   They  are  therefore  not  readily  transferable  to  a  different 
zonal  scheme  in  another  urban  area  unless  this  scheme  is  very  similar  in  size 
and  composition  to  the  initial  one. 

The  disadvantages  of  aggregate  zonal  regression  equations  are  of  even 
greater  concern  in  the  case  of  home-based  productions  (HBW  and  HBO) .   These 
equations  use  zonal  aggregate  data  and  assume  that  the  zones  are  homogeneous 
in  terms  of  the  socio-economic  characteristics  of  the  households  in  the  zone; 
in  many  instances,  this  might  not  really  be  the  case.  Aggregation  to  the  zonal 
level  tends  to  conceal  a  considerable  amount  of  the  variation  between  individual 
households  within  a  zone,  thus  masking  the  causal  nature  of  the  explanatory 
variables.   These  models  are  therefore  not  behavioral  in  nature,  and,  as  a 
consequence,  they  are  not  sensitive  to  changes  occurring  at  the  basic  decision 
level  of  tripmaking. 

Trip  Generation  Analysis  At  The  Disaggregate  Household  Level 
Within  the  framework  of  the  established  transportation  planning  process, 
two  types  of  disaggregate  household  trip  generation  models  can  be  used: 
regression  models  calibrated  with  household  data,  and  cross-classification 
(category  analysis)  models.   The  latter  stratify  households  according  to  their 
socio-economic  characteristics  and  provide  estimates  of  the  trip  rates  for  each 
of  the  household  categories.   These  models  have  numerous  advantages  over  zonal 
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aggregate  models.   In  addition,  category  analysis  offers  the  following  advantages 
over  the  regression  technique: 

-  ability  to  express  non-linear  relationships 

-  distribution-free  characteristic  whereby  no  limiting  assumptions  as  to 
the  statistical  distribution  of  the  variables  have  to  be  made  as  in  the 
case  of  regression  analysis 

-  ease  of  understanding  and  simplicity  of  application 

-  ease  of  monitoring  and  updating  using  small  disaggregate  samples  or  even 
subjective  judgment. 

These  reasons  make  cross-classification  models  preferable  to  regression  models, 
especially  for  simplified  transportation  planning  for  small  urban  areas. 

Cross-Classification  Procedure.   The  cross  classification  technique  for 
household  trip  production  analysis  stratifies  household  trip  rates  according  to 
the  major  household  socio-economic  characteristics  affecting  trip  making  behavior. 
Extensive  empirical  evidence  was  examined  to  indicate  that  the  two  major 
determinants  of  household  trip  generation  are  auto  ownership  and  household  size. 
Households  can  therefore  be  categorized  according  to  these  two  characteristics, 
and  average  trip  rates  for  each  household  category  (or  cell)  are  then  developed. 
In  order  to  do  so,  trip  rates  (by  trip  purpose)  for  the  individual  households 
are  computed  and  accumulated  within  each  cell.   The  average  number  of  trips 
made  by  households  as  well  as  the  standard  deviation  are  then  computed  for  each 
household  socio-economic  group.   These  rates  can  be  plotted  and  a  smooth  curve 
can  be  fitted  to  the  observations. 

Cross-Classified  Vehicle  Trip  Rates  for  Lafayette  and  Evansville.   Since 
appropriate  disaggregate  household  data  was  available  only  in  two  of  the  study 
areas,  Lafayette  and  Evansville,  household  cross-classified  trip  rates  were 
developed  only  for  these  two  urban  areas.   Trip  rates  were  determined  for  each 
purpose  (HBW,  HBO,  and  NHB) .   Also,  trips  were  accumulated  within  each  household 
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socio-economic  category,  and  average  trip  rates,  as  well  as  standard  deviations, 
were  computed  for  each  category.   In  Figures  1  and  2  are  shown  total  internal 
vehicle  trip  rates  for  Lafayette  and  Evansville,  respectively. 

Cell-by-cell  comparisons  of  trip  rates  by  purpose  between  the  two  urban 
areas  were  made  using  a  series  of  t-tests.   Only  a  few  cells  for  all  three 
purposes  exhibited  significant  differences  between  the  urban  areas.   Thus,  the 
comparisons  made  strengthened  the  belief  that  behavioral  disaggregate  household 
models  are  transferable  from  one  urban  area  to  another,  especially  in  the 
case  of  small  urban  areas,  regardless  of  the  socio-economic  differences  between 
them.   In  spite  of  several  such  differences  between  the  urban  areas  of  Lafayette 
and  Evansville,  the  household  trip  rates  were  not  significantly  different,  which 
means  that  the  trip  making  behavior  is  in  essence  similar  in  these  cities  even 
though  their  aggregate  travel  patterns  are  different. 

Bayesian  Updating  Of  Trip  Generation  Information 
Even  though  trip  rates  have  been  shown  to  be  stable  between  small  urban 
areas,  it  might  be  desirable  to  modify  them  so  that  they  reflect  local  conditions 
more  accurately.   This  can  be  done  by  updating  the  borrowed  parameters  using 
information  particular  to  the  area  under  study.   This  information  may  be  either 
objective  or  subjective. 

Bayesian  techniques  can  be  used  to  update  trip  rates  developed  by  category 
analysis  in  a  certain  area  before  applying  them  in  a  different  area.   In  this 
study,  the  methodology,  as  well  as  example  applications,  were  described  for  two 
cases: 

1.  New  information  obtained  from  small  disaggregate  sample  coupled  with 
the  prior  information  based  on  the  borrowed  travel  data  from  a  similar 
urban  area. 

2.  New  information  derived  from  subjective  judgment  coupled  with  the  prior 
information  based  on  the  borrowed  travel  data  from  a  similar  urban  area. 
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Conclusions  On  Trip  Generation  Phase 

I •   The  distribution  of  trips  by  purpose  is  significantly  different  between 
the  study  areas,  and  within  an  urban  area,  this  distribution  varies 
between  socio-economic  groups. 

2.  For  small  urban  areas,  differences  in  socio-economic  characteristics 
constitute  the  major  differentiating  characteristics  affecting  trip 
frequency;  the  socio-economic  mix  in  the  study  areas  was  indicated  by 
the  auto  ownership  distribution  between  households. 

3.  Caution  should  be  exercised  when  borrowing  zonal  trip  generation 
equations  for  use  in  urban  areas  other  than  the  one  where  they  were 
developed. 

4.  Home-Based  Work  zonal  production  equations  have  more  potential  for 
transferability  than  Home-Based  Other  equations  because  of  the  nature 
of  the  work  trip. 

5.  Non-Home  Based  zonal  production  equations  have  good  potential  for 
transferability  and  thus  can  be  used  successfully  in  synthetic 
applications. 

6.  Zonal  trip  attraction  regression  equations  can  be  used  to  synthesize 
zonal  attractions  in  areas  different  from  the  ones  where  they  were 
developed.  Non-Home  Based  and  Home-Based  Work  equations  predict  more 
accurately  than  Home-Based  Other  equations  even  though  the  latter 
predict  within  reasonable  limits.   However,  it  is  recommended  that 
the  synthesized  attractions  be  balanced  with  total  productions  as 
obtained  by  more  accurate  synthetic  techniques. 

7.  Disaggregate  household  trip  generation  rates  have  stronger  theoretical 
justification  for  transferability  because  they  are  behavioral  in  nature. 

8.  The  number  of  levels  needed  within  each  variable  in  the  cross-classification 
trip  generation  technique  has  been  determined  to  be: 
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-  Auto  ownership:   0  cars,  1  car,  2  cars,  3  or  more  cars. 

-  Household  size:   1,  2,  3,  4,  5  or  more  persons. 

9.   Bayesian  statistics  can  be  effectively  used  to  update  cross-classified 
household  trip  generation  rates.   Objective  or  subjective  information 
can  be  used  to  update  trip  rates  spatially  as  well  as  temporally. 

Recommended  Procedure  For  Synthetic  Trip  Generation 
Analysis  In  Small  Urban  Areas 
This  procedure  is  limited  to  internal  vehicle  trips.   Before  detailed  trip 
generation  analysis  is  undertaken,  rough  estimates  can  be  obtained  using  the 
trip  rates  given  in  Table  2  and  in  other  tables  in  Chapters  3  and  4  of  the 
report,  "Framework  for  the  Transferability  of  Trip  Generation  Parameters  of 
Small  Urban  Areas  in  Indiana"  (1) . 

The  following  steps  could  next  be  followed: 

a.  Establishing  Control  Totals 

This  could  be  done  using  the  disaggregate  trip  rates  as  shown  in  Figures  1 
and  2  whereby  areawide  control  totals  can  be  obtained  by  multiplying  the 
appropriate  rate  by  the  areawide  number  of  households  in  each  category.   More 
detailed  information  can  be  obtained  in  Reference  1. 

Control  totals  are  obtained  for  each  purpose  separately.   It  should  be 
noted  however  that  the  Non-Home  Based  trip  rates  reported  in  Chapter  6  of 
Reference  1  do  not  include  truck  trips.   The  missing  information  can  be  supplied 
either  by  using  the  areawide  vehicle  trip  rates  reported  in  Chapter  4,  or  by 
using  local  survey  information,  or  by  multiplying  the  number  obtained  from  the 
trip  rates  by  a  correction  factor  between  1.4  and  1.6  depending  on  the  area 
being  studied  (this  range  is  based  on  figures  obtained  from  the  study  areas) . 

b.  Trip  Productions 

Zonal  values  of  the  trip  productions  are  needed.   These  can  be  obtained  in 
the  following  way: 
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-  HBW  and  HBO:   using  appropriate  cross-classification  tables  and  multiplying 
by  the  zonal  number  of  households  of  each  category.   HBW  productions  can 
alternatively  be  determined  using  the  equation  recommended  for  that 
purpose  in  Chapter  5  of  Reference  1. 

-  NHB:   using  synthetic  equations  as  recommended  in  Chapter  5  of  Reference  1. 
The  above  zonal  productions  should  be  adjusted  to  agree  with  the  control 

totals  established  earlier. 

c.   Trip  Attractions 

Regression  equations  can  be  used  to  synthesize  zonal  attractions  provided 
that  the  total  attractions  are  balanced  by  total  productions  for  each  purpose. 
In  addition,  special  analyses  might  be  needed  for  special  generators  unique  to 
the  study  area. 

Alternatively,  attraction  rates  derived  from  national  data  can  be  used  as 
described  in  Chapter  5,  again  with  the  condition  that  the  total  attractions  are 
constrained  to  equal  the  total  productions. 

Finally,  the  whole  process  would  have  to  be  "tuned"  so  that  the  volumes 
assigned  to  the  network  (after  the  trip  distribution  and  trip  assignment  steps 
are  executed)  compare  favorably  with  the  observed  volume  counts. 

SYNTHETIC  TRIP  DISTRIBUTION 
The  second  phase  of  the  study  concentrated  on  the  trip  distribution  part 
of  the  demand  modelling  process  (2) .   The  performance  of  some  of  the  various 
suggested  synthetic  gravity  model  calibration  techniques  were  evaluated  and  an 
improved  procedure  for  synthetic  self-calibration  of  the  gravity  model  was 
developed.   Data  from  a  wide  range  of  small  urban  areas  were  used  in  determining 
variability  of  different  travel  characteristics.   In  addition,  detailed  travel 
data  from  three  urban  areas  in  Indiana,  Lafayette,  Anderson  and  Muncie,  were 
used  in  the  development  of  the  recommended  procedure. 
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Travel  Time  Characteristics 

In  Table  3  are  given  the  relevant  travel  time  characteristics  for  each  of 
the  study  areas.  A  detailed  analysis  was  conducted  to  evaluate  the  travel 
time  characteristics  and  it  was  concluded  that  the  prediction  of  average 
travel  time  based  on  the  population  size  of  the  urban  area  alone  would  be  of 
little  practical  use.   This  conclusion  was  derived  from  the  fact  that  even 
though  a  positive  correlation  existed  between  average  travel  time  and  popula- 
tion size,  the  variance  of  the  data  about  the  regression  line  was  substantial 
as  measured  by  low  R-squared  values.   This  results  in  a  regression  equation 
with  a  confidence  interval  too  large  to  be  of  practical  use. 

A  comparison  of  home  based  work  travel  time  frequency  distributions  for 
the  Indiana  study  areas  suggested  that  more  careful  selection  of  zonal  terminal 
times  might  enable  better  comparisons  to  be  made  between  the  various  study  areas. 
Even  though  two  distributions  may  look  similar,  however,  the  average  travel 
times  can  differ  significantly.   Using  travel  time  frequency  distributions  of 
one  urban  area  to  compare  with  the  distributions  generated  from  another  urban 
area's  synthetic  trip  distribution  model  is  therefore  of  little  benefit,  unless 
it  is  previously  known  that  the  average  travel  time  in  the  two  areas  are  very 
close. 

The  above  points  verify  the  important  fact  that  travel  time  statistics  of 
an  urban  area  are  a  unique  result  of  the  area's  characteristics  and  are  there- 
fore difficult  to  compare  with,  or  use  for  prediction  in  another  study  area. 

Gravity  Model  Sensitivity  Analysis 
Standard  Gravity  Model  Formulation.   The  gravity  model  has  been  the  most 
widely  used  trip  distribution  technique,  mainly  because  it  is  simple  in  concept 
and  the  procedure  is  well  documented.  Mathematically,  the  gravity  model  is 
stated  as  follows: 
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where: 


T   =  trips  produced  in  zone  i  and  attracted  to  zone  j, 

P.   =  trips  produced  by  zone  i, 

A.   =  trips  attracted  by  zone  j,  and 

F  .  =  friction  factor  which  expresses  the  effect  of  spatial  separation 

on  trip  interchanges. 
K..  =  zone-to-zone  adjustment  factors.   (These  factors  are  ignored  in 

synthetic  techniques  because  origin-destination  data  is  needed  to 

determine  these  values.) 
Friction  factors  can  be  represented  by  a  set  of  values  at  one  minute  time 
increments  or  by  an  equation  of  the  form: 

F(t)  =  Ctae6t 
where: 

F(t)  =  friction  factor  at  time  t, 

C    =  constant  of  proportionality, 

a         =  alpha,  calibration  parameter, 

3    =  beta,  calibration  parameter,  and 

e    =  base  of  the  natural  logarithms. 
Of  particular  interest  to  this  study  is  the  alternate  form  of  this  equation 
known  as  the  negative  exponential  function  where: 

F(t)  =  Ce"6t 
with  the  variables  defined  as  above. 

Gravity  Model  Calibrations.   The  calibration  of  a  gravity  model  basically 
involves  repeated  adjustments  of  a  set  of  friction  factors  until  the  generated 
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gravity  model  travel  time  frequency  distribution  replicates  the  0-D  survey 
travel  time  distribution  within  a  predetermined  error  limit.   This  procedure, 
however,  requires  that  an  0-D  survey  has  been  performed  in  order  to  conduct 
the  accuracy  tests  and  the  calibration  itself. 

In  order  to  eliminate  the  cost  and  time  involved  in  conducting  an  0-D 
survey,  researchers  and  transportation  planners  have  occasionally  borrowed 
friction  factors  from  other  urban  areas  possessing  similar  characteristics,  and 
then  applied  them  to  a  new  study  area.  An  analysis  was  performed  to  evaluate 
the  performance  of  the  gravity  model  utilizing  various  synthetic  techniques  to 
assess  possible  errors  that  might  be  introduced  by  the  borrowing  of  friction 
factors  from  a  similar  urban  area  or  by  the  use  of  standardized  friction  factors 
or  other  synthetic  techniques.   The  transportation  study  areas  of  Lafayette  and 
Anderson  were  used  for  the  purpose  of  this  evaluation.   The  population  sizes 
of  these  two  areas  are  approximately  equal,  but  spatial  activity  distributions 
and  network  impedances  are  quite  different  from  each  other. 

This  portion  of  the  research  consisted  of  four  parts.   First,  tests  were 
conducted  to  measure  the  sensitivity  of  the  gravity  model  with  respect  to  small 
errors  in  the  friction  factors.   Simultaneously,  these  tests  revealed  the 
sensitivity  of  the  gravity  model  with  respect  to  the  network  characteristics 
and  land  use  activity  distributions  of  the  individual  study  areas.   Second, 
the  variation  of  friction  factors  with  respect  to  population  size  was  investigated, 
This  investigation  involved  the  examination  of  the  concept  of  similar  urban  areas 
from  which  friction  factors  may  be  borrowed.   Third,  an  evaluation  was  performed 
using  standardized  a   parameters  in  the  gravity  models  for  the  Lafayette  and 
Anderson  urban  areas.   Fourth,  the  friction  factors  of  Lafayette  and  Anderson 
were  exchanged  or  borrowed  from  each  other  in  order  to  test  how  borrowed  friction 
factors  may  affect  trip  distribution  modelling  accuracies. 
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Sensitivity  analysis  of  the  gravity  models  for  Lafayette  and  Anderson 
resulted  in  three  important  observations  in  regards  to  synthetic  trip  distribu- 
tion modelling. 

1.  Acceptable  model  results  will  occur  when  the  borrowed,  or  standardized 
friction  factors  are  very  similar  to  those  that  would  have  been 
normally  calibrated  from  0-D  survey  data. 

2.  Acceptable  model  results  may  occur  when  synthetic  methods  are  used 
because  the  urban  area's  gravity  model  is  not  sensitive  to  friction 
factor  errors. 

3.  Acceptable  model  results  may  occur  when  one  individual  trip  purpose 
error  is  offset  by  an  opposite  error  in  one  or  more  of  the  other  trip 
purposes. 

Clearly,  the  first  situation  is  the  most  preferred  with  the  second  also 
being  acceptable.   The  third  situation,  however,  is  unacceptable  because 
chance  plays  a  larger  role  in  the  model's  outcome  than  the  procedure  by  which 
the  friction  factors  were  obtained.   In  addition,  analysis  of  data  from  a 
number  of  small  urban  area  transportation  studies  resulted  in  one  other 
important  conclusion:   population  size  is  an  inadequate  criterion  by  which  to 
select  a  similar  urban  area  for  borrowing  friction  factors  or  their  parameters. 

Development  and  Testing  Of  Procedures  For  Synthetic 
Self-Calibration  Of  Gravity  Models 
Areawide  Approach  to  Synthetic  Self-Calibration  of  a  Gravity  Model.   It  is 
possible  to  relate  mathematically  the  parameter  of  a  friction  factor  equation  to 
the  spatial  land  use  activity  distribution  and  network  impedances  through  the 
use  of  the  areawide  descriptors  of  mean  travel  time  and  travel  time  variance. 
This  is  because  the  friction  factor  equation  of  a  gravity  model  represents  the 
trip  makers'  destination  choice  process  with  respect  to  the  spatial  distribution 
of  the  land  use  activities  as  represented  by  the  destination  choice  variable  of 
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interzonal  travel  time.   The  aggregate  effect  of  all  the  trip  makers'  destina- 
tion choices  result  in  a  total  travel  pattern  that  can  be  approximately 
described  in  terms  of  the  areawide  travel  time  statistics  of  mean  and  variance. 

Data  used  in  the  study  of  the  areawide  approach  consisted  of  the  8 
parameters  of  the  negative  exponential  friction  factor  function,  areawide 
average  travel  times  and  travel  time  variances,  study  area  population,  and 
study  area  density  as  computed  or  reported  for  eighteen  small  urban  area 
transportation  studies.   Regression  analysis  was  used  to  search  for  mathematical 
relationships  between  the  8  parameters  and  the  areawide  travel  time  statistics 
and  any  other  variables  of  possible  interest. 

The  self -calibrating  8  parameter  procedure  using  the  regressions  equations 
was  found  to  be  not  as  reliable  a  procedure  as  the  traditionally  calibrated 
gravity  model  technique.   Factors  related  to  the  poor  performance  of  this 
procedure  were  considered  in  the  development  of  the  second  approach  to  synthetic 
self -calibration  of  a  gravity  model. 

Origin  Zone  Specific  Approach  to  Synthetic  Self-Calibration  of  a  Gravity 
Model.   The  proposed  synthetic  procedure  involves  an  origin  zone  specific  self- 
calibrating  model  in  which  the  only  input  data  required  are  the  zonal  productions 
and  attractions,  a  zone-to-zone  travel  time  matrix  (skim  tree),  and  the  origin 
zone  terminal  times.   A  travel  time  distribution  is  no  longer  needed  for 
calibration,  thereby  eliminating  the  necessity  for  an  0-D  survey  for  internal- 
internal  trips. 

The  study  areas  of  Lafayette,  Anderson,  and  Muncie,  Indiana,  were  used  to 
generate  a  data  base  because  of  data  availability  and  the  wide  sampling  space 
provided  by  the  differing  characteristics  of  these  areas.   The  gravity  models 
of  these  three  areas  were  calibrated  at  the  zonal  level  to  generate  the  following 
information  for  each  origin  zone  and  for  each  trip  purpose  -  HBW,  HBO,  and  NHB 
trips: 
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8.    =  negative  exponential  friction  factor  parameter  for  zone  i  of  the 

origin  zone  calibrated  gravity  model, 
0.    =  trips  originating  from  zone  i, 
AVEO.  =  opportunity  average  travel  time  for  origin  zone  i  when  8  equals 

0.0,  that  is,  friction  factor  of  constant  value  1.0, 
VARO.  =  opportunity  travel  time  variance  for  origin  i  when  8.  equals  0.0, 
AVE.   =  origin  specific  average  travel  time  for  zone  i  of  calibrated 

gravity  model  (same  as  0-D  survey  value) , 
VAR.   =  origin  specific  travel  time  variance  for  zone  i  of  calibrated 

gravity  model,  and 

TERM.  =  origin  zone  terminal  time. 

The  above  variables  were  then  used  in  a  curve  fitting  regression  procedure  so  as 

to  find  a  relationship  of  the  following  form: 

8 .  =  f (AVEO . ,  AVE . ,  VARO . ,  VAR . ,  TERM . ) . 
l         l     l      i     i      l 

The  initial  regression  analysis  resulted  in  an  equation  form  possessing 

very  acceptable  statistics  for  the  curve  fitting  of  the  origin  zone  specific 

8  parameters  to  their  respective  zonal  travel  time  statistics  of  mean  and 

variance.  However,  inspection  of  the  residual  plot  revealed  a  greater  than 

expected  number  of  outliers.   Careful  review  of  the  zonal  data  and  a  sensitivity 

analysis  led  to  the  following  modifications: 

-  removal  of  the  variable  VAR.  from  the  analysis,  and 

-  elimination  of  the  simultaneous  inclusion  of  both  AVE.  and  AVEO.^  in  the 
regression  equation. 

A  high  correlation  observed  between  AVEO.  and  AVE.  indicated  that  an  indirect 
technique  to  estimating  the  8  parameter  could  be  possible.   Consequently,  a 
revised  procedure  was  developed  with  AVEO.  as  the  key  variable. 


25 


Regression  analysis  of  AVE  on  AVEO  and  TERM  was  accomplished  by  weighting 

the  zonal  data  in  order  to  properly  balance  the  data  by  trip  productions  (origins) 

and  by  trip  purpose.   The  regression  analysis  result  for  trip  purpose  HBW 

follows: 

AVE,,   =   -0.4328  +  1.2761   TERM.   +  0.8921    (AVEO.    -  TERM.) 
i  ill 

R2  =    .865 

Similar  equations  for  HBO  and  NHB  trips  were  found,  and  all  equations  possessed 

R-squared  values  higher  than  those  achieved  in  earlier  attempts.   The  strong 

statistical  fit  of  the  equations  shows  the  importance  of  AVEO.,  the  opportunity 

average  travel  time,  on  the  destination  choice  process. 

Incorporation  of  the  Estimating  Equations  into  a  Gravity  Model.   The 

recommended  procedure  for  synthetic  calibration  was  based  on  three  assumptions 

regarding  the  interaction  of  origin  zone  average  travel  time  and  the  zonal 

3  parameter: 

-  The  origin  zone  average  travel  time  varies  linearly  with  the  origin  zone 
3  parameter. 

-  Zonal  average  travel  time  is  independent  of  the  B  parameters  of  other 
zones. 

-  Zonal  average  travel  time  is  independent  of  the  attraction  constraint 
iteration  procedure. 

The  estimating  equations  of  AVE.,  zonal  average  travel  time,  and  the  three 
assumptions  listed  above  were  incorporated  into  a  gravity  model  procedure 
described  below  in  conjunction  with  Figure  3. 
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STEP  I 


0.0 


Step  1  -  Average  at  ft  equal  zero 
(0.0)  -  AVEO 

Step  2  -  Average  at  ft  equal  to  a 
set  value. 

Step  3  -  Average  from  regression 
equation,  AVE 

Step  4  -  Interpolate  for  new  ft 


P> 


SET  VALUE 


Origin  Zone  ft  Parameter 


Figure  3.   Generalized  Flow  Chart  of  Synthetically  Self-Calibrated, 

Origin  Zone  Specific  Gravity  Model  as  Applied  to  a  Typical 
Zone 
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First,  trips  are  distributed  with  all  zonal  6  parameters  equal  to  0.0, 
and  then  the  resulting  zonal  opportunity  average  travel  times,  AVEO  ,  are 
computed.   Second,  trips  are  distributed  with  all  origin  zone  6  parameters 
equal  to  a  preselected  value,  and  then  the  origin  zone  average  travel  times 
are  computed.   Third,  using  the  estimating  equation  for  the  appropriate  trip 
purpose  and  the  AVEO.  values  computed  in  step  one,  the  origin  zone  average 
travel  times,  AVE  ,  are  determined.   Fourth,  using  the  data  points  from  steps 
one  and  three  and  the  estimated  origin  zone  average  travel  times,  AVE.,  the 
zonal  3  parameters  are  obtained  through  interpolation.   All  trips  are  then 
distributed  using  the  interpolated  origin  zone  specific  8  parameters. 

Tests  conducted  on  three  separate  study  areas  indicated  that  the  proposed 
model  is  able  to  reproduce  trip  patterns  as  accurately  as  the  traditionally 
calibrated  gravity  model  procedure  which  is  calibrated  on  0-D  survey  data. 
The  opportunity  average  travel  time,  AVEO.,  is  the  key  variable  on  which  the 
model  was  developed  and  upon  which  the  very  acceptable  model  results  are  due. 

Conclusions  On  Trip  Distribution  Phase 

1.  The  need  and  importance  of  friction  factors  was  established  by  the 
values  of  average  travel  time  and  traffic  assignment  statistics  gen- 
erated by  trip  distribution  models  without  such  time  dependent  factors. 

2.  The  sensitivity  of  the  gravity  model  to  friction  factor  errors  may 
be  highly  dependent  upon  the  particular  urban  area  being  modelled. 

3.  A  large  variability  exists  in  the  correlation  between  population  size 
and  friction  factors. 

4.  Because  the  variability  of  average  travel  time  with  respect  to  popula- 
tion size  is  large,  the  validation  of  a  synthetic  trip  distribution 
model  by  comparing  only  the  average  travel  times  is  not  a  valid 
approach. 
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5.  While  general  shapes  of  travel  time  frequency  distributions  between 
two  areas  may  appear  quite  similar,  even  a  small  variation  or  shift 
in  a  distribution  can  significantly  alter  average  travel  time  and 
therefore  also  the  accuracy  of  the  model. 

6.  The  most  reliable  synthetic  trip  distribution  modelling  technique  is 
the  origin  zone  specific  synthetically  self-calibrating  gravity  model. 

7.  The  most  important  feature  of  the  proposed  procedure  is  its  more 
direct  linkage  to  the  spatial  distribution  of  land  use  activities 
within  an  urban  area  as  quantified  by  the  origin  zone  specific 
opportunity  average  travel  times. 

8.  No  collection  of  0-D  survey  data  is  needed  for  this  procedure,  nor 
is  it  necessary  to  find  a  similar  urban  area  from  which  to  borrow 
information  or  to  which  model  validation  comparisons  are  to  be  made. 

Computer  Program  For  The  Recommended  Procedure 
A  computer  program  JMGM  was  written  for  the  purpose  of  using  the 
recommended  procedure  with  the  UMODEL  program  of  the  UTPS  package.   A  PLANPAC 
version  of  JMGM  was  also  coded.   A  FORTRAN  IV  program  is  also  available  for 
incorporation  into  other  packages.   In  Appendix  is  given  a  listing  of  this 
program. 

EFFECTS  OF  SIMPLIFYING  ZONE  AND  NETWORK  SYSTEMS  ON  THE 
ACCURACY  OF  TRAFFIC  ASSIGNMENT 

The  third  phase  of  the  study  included  an  evaluation  of  the  effects  of 

simplifying  zone  and  network  systems  on  the  accuracy  of  traffic  assignment  (3). 

major  factors  of  interest  in  this  phase  were  average  traffic  zone  size,  use  of 

census  tracts  as  traffic  zones,  level  of  detail  used  in  street  networks,  and 

method  of  traffic  assignment  used.   Since  a  statistical  design  of  experimentation 

was  needed  to  test  these  factors  and  the  possible  interactions  among  them, 
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actual  transportation  network  and  travel  data  from  two  small  urban  areas  in 
Indiana,  Lafayette  and  Anderson,  were  used. 

Alternative  Traffic  Zone  Systems.   For  each  urban  area,  alternative  traffic 
zone  systems  were  developed  by  substituting  existing  zones  with  smaller  numbers 
of  large  zones.   Only  the  internal  zones  were  changed;  all  external  zones  were 
left  intact.   Larger  zones  were  developed  by  combining  existing  zones  into 
larger  units.   Only  adjacent  existing  zones  with  similar  land  uses  were 
combined  into  a  larger  zone.   The  numbers  of  vehicle  trip  productions  and 
attractions  (P's  and  A's)  for  a  new  zone  were  computed  simply  by  adding  to- 
gether the  P's  and  A's  of  its  component  zones.   New  zone  systems  were  also 
developed  from  census  tracts. 

When  a  new  zone  was  created,  the  centroid-connecting  links  to  the  old 
zones  were  changed.   Only  centroid  connectors  that  were  near  the  new  zone 
centroid  and  collectively  provide  access  to  all  nearby  arterial  and  collector 
streets  were  chosen  for  the  new  zone. 

Terminal  Time  and  Intrazonal  Travel  Time  Estimations.   New  terminal  times 
and  intrazonal  travel  times  were  established  for  each  new  zone.   The  terminal 
time  for  a  new  zone  was  set  as  the  terminal  time  predominant  among  the  existing 
zones  from  which  the  new  zone  was  formed.   The  following  equation  was  developed 
to  compute  the  intrazonal  travel  times  for  new  zones: 
I  (ITT  *  p  ) 

ITT-w  =  TF=-  z=1>2 ■  n 

Jz  pz 

where: 

ITT    =  intrazonal  travel  time  for  new  zone, 
new 

ITT    =  intrazonal  travel  time  for  old  zone, 

n      =  no.  of  old  zones  within  new  zone,  and 

P      =  areal  function  of  old  zone  Z  within  new  zone. 
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Alternative  Network  Systems.   Two  alternative  street  network  syste 
were  developed  for  each  urban  area,  one  being  the  full,  existing  coded  network 
for  an  urban  area.   The  second  alternative  was  a  reduced  version  of  the  full 
network,  with  all  links  on  collector  streets  that  carried  less  than  5000 
vehicles  per  day  removed.   All  links  on  arterial  streets  and  high-volume 
collector  streets  were  kept  in  the  reduced  networks,  as  well  as  in  the  full 
networks,  since  this  level  of  street  network  detail  was  a  minimum  level  at 
which  all  important  links  would  still  be  represented. 

Procedures  For  Data  Collection  And  Analysis 
For  both  Lafayette  and  Anderson,  each  alternative  traffic-zone  system 
was  tested  with  each  alternative  street  network.   This  resulted  in  ten  zone- 
networks  tests  for  Lafayette  and  eight  for  Anderson.  Each  test  included  a 
gravity  model  calibration,  an  all-or-nothing  assignment,  and  a  capacity  re- 
straint assignment.   After  all  tests  had  been  completed,  the  data  resulting 
from  the  calibrations  and  assignments  were  analyzed. 

In  conducting  the  zone-network  tests  considerable  street  network  and 
travel  survey  data  were  required.   The  street  network  data  included  a  listing 
of  all  links  with  each  link's  length,  travel  time,  ground  count,  and  capacity. 
Travel  survey  data  included  vehicle-trip  P's  and.  A's  for  each  zone,  trip 
length  distributions,  and  vehicle  trips  interchanged  between  external- internal 
.and  external-external  aone. pairs. 

Each  test  followed  the  same  set  of  procedures.   These  activities  were: 

-  Preparation  of  network  and  travel  data. 

-  Gravity  Model  calibration. 

-  Traffic  assignment. 

The  final  activity  involved  as  its  last  step  the  comparison  of  assigned  volumes 
with  ground  count  volumes.   This  step  produced  the  data  on  which  most  of  the 
analyses  were  based. 
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Analyses  of  Data.   The  data  from  the  test  calibrations  and  assignments 
were  analyzed  by  several  methods.   Aggregate  methods  were  used  to  evaluate 
overall  model  calibrations  and  traffic  assignments.   Stratified  methods  were 
used  to  study  traffic  assignments  in  more  detail. 

Alpha  values  and  intrazonal  trip  measurements  were  used  in  the  aggregate 
analysis  of  model  calibrations  and  trip  distributions  of  all  tests.   Alpha 
values  were  examined  for  their  stability  with  changes  in  the  zone  system  or 
street  network.   Intrazonal  trips,  measured  as  a  percentage  of  all  internal 
trips,  could  be  expected  to  increase  as  zone  size  increased.   Total  vehicle- 
miles  travelled  (VMT)  and  overall  percent  root  mean-square  error  (%  RMSE) 
were  aggregate  measures  for  comparing  traffic  assignments. 

To  study  traffic  assignments  at  a  detailed  level,  the  total  population 
of  links  in  a  street  network  was  stratified  into  volume  groups  based  on  ground 
count  volumes . 

Statistical  analysis  at  the  stratified  level  was  conducted  in  two  phases. 
In  the  first  phase,  the  %  RMSE's  of  volume  groups  in  each  assignment  were 
analyzed  for  the  main  effects  and  certain  interactions.   The  interests  lay 
with  the  main  effects  of  and  possible  interactions  between  certain  factors: 
the  number  of  traffic  zones,  the  type  of  zone  boundary,  the  type  of  street 
network,  the  method  of  traffic  assignment,  and  the  volume  group.   Two-way 
interactions  were  examined  and  multiple-factor  analyses  of  variance  (ANOVA's) 
were  run  and  compared  to  determine  main  effects  and  interactions. 

For  the  second  phase  of  stratified  analysis,  observations  were  the 
absolute  differences  between  the  assigned  volumes  and  ground  counts  on 
individual  links.   Links  on  a  network  are  not  independent  of  each  other  and 
neither  are  their  traffic  volumes.   Without  independent  observations,  para- 
metric statistical  methods  lack  validity.   For  this  reason,  these  data  were 
analyzed  by  nonparametric  methods  using  Freidman  rank  sums. 
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Analyses  Of  Test  Results 

Definite  patterns  appeared  as  to  how  zone  and  network  factors  affect 
values  of  the  alpha  parameter  in  gravity  model  calibration.   Reducing  the 
number  of  zones  seemed  to  reduce  the  alphas  for  HBW  trips,  while  the  alphas 
of  the  other  trip  purposes  were  only  slightly  affected.   Reducing  the  number 
of  links  in  the  network  seemed  to  increase  alphas  for  all  trip  purposes. 
Alpha  values  were  more  variable  for  Lafayette  than  for  Anderson.   This 
occurrence  may  have  been  due  to  differences  in  the  characteristics  of  the 
two  cities,  to  differences  in  the  accuracy  of  the  data  from  the  two  cities, 
or  to  the  use  of  very  large  zones  in  Lafayette. 

The  effects  of  zone  and  network  factors  on  overall  VMT  appeared  small. 
Total  VMT  in  capacity-restraint  assignments  generally  varied  less  than  5% 
with  changes  in  zone  size.   Reducing  network  detail  increased  total  VMT  for 
both  cities. 

Decreasing  the  number  of  zones  tended  to  increase  the  %  RMSE.   Assigning 
traffic  by  capacity-restraint  always  reduced  the  %  RMSE.   The  use  of  Capacity- 
restraint  reduced  the  effect  of  the  number  of  zones  on  the  %  RMSE.   Reduced 
networks  generally  showed  lower  %  RMSE's,  while  no  real  trend  could  be 
established  for  the  effect  of  zone  boundary  type  on  overall  %  RMSE. 

In  the  stratified  analysis  of  %  RMSE,  three  ANOVA's  were  run  at  a  .05 
level  of  significance.   Observations  were  the  %  RMSE's  by  volume  group  for 
every  traffic  assignment.   The  ANOVA  showed  that  three  factors,  number  of 
zones,  assignment  method,  and  volume  group,  had  particularly  significant  main 
effects  on  the  accuracy  of  traffic  assignments.   Interactions  between  factors 
were  few  and  generally  inconsistent.   The  nonparametric,  link-by-link  analyses 
describe  the  main  effects  more  fully. 
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The  results  of  these  analyses  implied  that  several  factors  have  effects 
on  the  accuracy  of  traffic  assignments: 

-  Using  fewer  (and  larger)  zones  tends  to  reduce  the  accuracy  of 
traffic  assignments; 

-  Capacity-restraint  assignment  is  substantially  more  accurate  than 
all-or-nothing  assignment; 

-  Limited  reduction  in  the  detail  of  a  street  network  will  not  materially 
affect  the  accuracy  of  traffic  assignments  on  high  volume  streets; 

-  The  boundaries  of  census  tracts  will  serve  as  a  basis  for  boundaries 
of  traffic  zones  without  adverse  effects  on  assignment  accuracy. 

Guidelines  For  Zone  System  And  Network  Development 
The  purpose  of  this  research  was  to  study  the  effects  of  reducing  the 
number  of  traffic  zones  and  simplifying  street  networks  on  traffic  assignment 
accuracy.   The  following  proposed  guidelines  for  establishing  zone  systems  and 
network  configurations  in  emerging  metropolitan  areas  are  the  final  results  of 
this  endeavor. 

Zone  size  is  a  significant  factor  in  the  accuracy  of  traffic  assignments. 
The  average  zone  size  for  a  transportation  study  area  should  not  exceed  3  sq. 
mi.,  and  no  zone  should  exceed  15  sq.  mi.  in  area.   Otherwise,  gross  instability 
in  friction  factor  parameters  may  occur.   Use  of  90  internal  zones  is  suggested 
for  small  urban  areas  such  as  Lafayette  and  Anderson  in  order  to  achieve 
acceptable  overall  accuracy  levels. 

Census  tract  boundaries  can  be  used  successfully  as  a  basis  for  establishing 
a  traffic  zone  system.   The  inconsistencies  that  may  occur  between  census  tract 
boundaries  and  land  use  boundaries  have  little  effect  on  traffic  assignment 
accuracy. 
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A  reduced  street  network,  comprised  only  of  collectors  with  an  ADT  over 
5000  and  all  higher-class  streets  and  highways,  should  be  adequate  for 
planning  purposes. 

Finally,  capacity-restraint  traffic  assignment  is  suggested  for  assigning 
trips  to  accurately  simulate  ground  counts.   Not  only  is  capacity-restraint 
assignment  more  accurate  than  all-or-nothing  assignment,  but  capacity-restraint 
assignment  also  reduces  the  adverse  impact  that  the  use  of  a  small  number  of 
zones  can  have  on  assignment  accuracy. 

CONCLUSIONS 

In  this  report  was  summarized  the  findings  of  the  project  on  synthetic 
demand  modeling  techniques  in  transportation  planning  for  urban  areas  in  Indiana. 
Only  internal  vehicle  trips  were  considered  and  guidelines  were  developed  for 
trip  generation,  trip  distribution  and  traffic  assignment  parts  of  travel  demand 
forecasting. 

Efforts  are  being  made  by  the  planning  staff  of  the  Indiana  State  Highway 
Commission  to  implement  the  study  results  in  several  urban  areas  in  Indiana. 
At  present  implementation  in  three  urban  areas  are  being  undertaken. 

Some  of  the  results  of  this  project  are  being  used  in  Bloomington, 
Indiana.   In  particular,  the  internal  trip  generation  information  for  Bloom- 
ington is  being  developed  on  the  basis  of  the  study  results.   Plans  are  being 
made  to  use  the  trip  distribution  part  in  the  Lake-Porter  (NIRPC)  area  as 
soon  as  a  network  becomes  available.   Use  of  the  guidelines  in  the  South 
Bend  (MACOG)  area  will  include  both  the  distribution  and  assignment  phases. 

It  is  felt  that  the  procedures  developed  in  the  study  will  greatly 
reduce  the  cost  and  time  involved  in  data  collection  and  model  development 
efforts  associated  with  conventional  travel  demand  modeling  process  based  on 
origin-destination  survey.   The  information  can  be  used  either  in  connection 
with  a  new  transportation  study  for  an  emerging  metropolitan  area  or  for 
updating  of  existing  transportation  plans. 
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APPENDIX 
COMPUTER  PROGRAM  FOR  SYNTHETIC  CALIBRATION  OF  GRAVITY  MODEL 
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THIS  IS  THE  CODE  FOR  ENTRY  POIHT  M0D13B 


ZONAL  UARIABLE3. 


(OS 
(OR 


PROGRAMMED  BY  JIM  MEKEMSOM.  PURDUE  UNIUEPSITY,  SEPT.  1373         JMGM 
IMPLEMENTED  AT  THE  INDIANA  STATE  HIGHWAY  COMMISSION  BY  NIKE  FARPISJMGM 

JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
C  JMGM 


c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C 

c 

C 

c 
c 
c 
c 

c 
c 
c 

c 
c 
c 


1 

2 

3 

4 

5 

S 

7 

8 

9 

10 

11 

12 

13 

14 


ZCNE  NUMBER 

ZONAL  PRODUCTIONS 

ZONAL  ATTRACTIONS 

TERMINAL  TIME.  ORIGIN  ZONE 

INTRAZONAL  TRAUEL  TIME 

ZONAL  BETA(IZ) 

GRAUITY  MODEL  DENOMINATORS 

NEW  ITERATED  UALUE  OF  Af J) 

SUM  COMPUTED  A(J)  =  SUM(IZ)  T(I,J) 

ZONAL  AUE.  AT  BETA=0.0 

ZONAL  UAR.  AT  BETA=0.0 

ZONAL  AUE.  AT  BETA=PBETA(TP)  OR  FINAL  UALUE 

ZONAL  UAR.  AT  BETA=PBETA(TP)  OR  FINAL  UALUE 

USED  ONLY  AS  WORKING  STORAGE.   TEMPORARY  FRICITON  FACTOR 


ORIGINS) 
DESTINATIONS) 


:) 


INITIALIZE  ZONAL  MATRIX 

IF((ITN0.NE.1).0R.(IZ.NE.1))RETURN 
DO  30  I=5rl3 
DO  30  J=l,OMAX 
30  Z(I.J)=0.0 


THIS  IS  THE  CODE  FOR  ENTRY  POINT  M0D13C 


REAL«4  PEETA(3).FCTR(99) 
PBETA(1)=0.15 
PBETA(2)=0.30 
PBETA(3)=0.30 


»«J"GM 

JMGM 
JMJM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 


C 
C 

c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


THIS  IS  THE  CODE  FOR  ENTRY  POINT  MOD13D 

ORIGIN  ZONE  SPECIFIC  SYNTHETICALLY  SELF-CALIBRATING  GRAUITY  MODEL 

INTERNAL-INTERNAL  TRIPS  ONLY4- 


USER  PARAMETER  INPUT  UIA  UPARMS 

UPARMS(l)  =  1   HBW 

2  HBO 

3  NHB 

UPARMS (2)  =  1   COMBINED  REGRESSION  EQUATIONS 
2   LAFAYETTE  REGRESSION  EQUATIONS 


JMGM 
JMGM 
JMCM 
JVGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 


1 

2 

1 
4 
5 
6 

7 

a 

3 
10 

li 

12 
13 

14 
15 
16 
17 
13 
19 
20 
21 
22 
23 
24 
25 
29 
27 
28 
29 
30 
31 
32 
33 
34 
35 
3S 
37 
33 
39 
40 
41 
42 
43 


46 
47 
43 
49 
50 
51 
52 
53 
54 
55 
55 
57 
5S 
5S 
GO 
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C 
C 
C 

c 
c 

c 

c 
c 


ANDERSON  REGRESSION  EQUATIONS 
NUNCIE  REGRESSION  EQUATIONS 


C 

c 

c 


c 
c 
c 


101 
20 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


21 


100 
50 


99 


c 

c. 

c. 

c 

c 

c 


NOTE:   OMAX=DMAX=  NUMBER  OF  INTERNAL  ZONES  ASSUMED. 

REAL»4  AUE.  UAR, ATT.  TUAR. SUM,  TRIP, RAUE. SX, SXX, TOT 

INITIALIZE  OR  CONTINUE 

IF((ITN0.EQ.1).AND.(IZ.EQ.1))G0  TO  101 
GO  TO  100 

COMPUTE  TOTAL  TRIPS 

TOT=0.0 

DO  20  J=l,OMAX 

T0T=T0T+Z(2.J) 

INITIALIZE 

SX=0.0 

SXX=0.0 

ATT=0.0 

TUAR=0.0 

CITER=0 

AITER=2 

DO  21  J=l,OMAX 

Z(8.J)=Z(3.J) 

Z(7.J)=T0T 

ITERATION  FLAGS  FOR  DECISION  TREE. 

CITER=CALIBRATION  FLAG 

0  -  BETA=0.0 

1  -  BETA=PBETA(TP) 

2  -  BETA=CALIBRATED  UALUE  FROM  REGRESSION  EQUATIONS 

AND  INTERPOLATION. 

AITER=ATTRACTION  FLAG 

0  -  BEGIN  ATTRACTION  ITERATION 

1  -  CONTINUE  ATTRACTION  ITERATION 

2  -  ATTRACTION  ITERATION  COMPLETE,  I.E.  CLOSURE 

IFCIZ.EQ.1)URITE(S,50)ITNO.CITER,AITER 

FORMAT(10X,7HITNO  =  , 15, 10X.8HCITER  =  , 15, 10X.8HAITER  =  ,15) 

IF((CITER.EQ.l).AND.CAITER.EQ.0).AND.(IZ.EQ.l))GO  TO  200 

IF((CITER.EQ.2).AND.(AITER.EO.0).AND.(IZ.EQ.l))GO  TO  300 

IF(AITER.£Q.0)GO  TO  400 

IF(AITER.EQ.l)GO  TO  500 

GO  TO  S00 


CITER=1,AITER=0,IZ=1    SET  BETACD'S  AND  FRICTION  FACTORS  TO 

200  TP=UPARMS(1) 
DO  201  J=l,99 

201  FCTR(J)=EXP(-PBETA(TP)»J) 
GO  TO  99 


JMGM  El 

JMGM  62 

JMGM  G2 

JMGM  G4 

JMGM  G5 

JMGM  55 

JMGM  57 

JMGM  58 

JMGM  E9 

JMGM  70 

jmgm  71 

JMGM  72 
JMGM  73 
JMGM  74 
JMGM  75 
JMGM  75 
JMGM  77 
JMGM  73 
JMGM  79 
JMGM  80 
JMGM  81 
JMGM  82 
JMGM  83 
JMGM  £4 
JMGM  35 
JMGM  35 
JMGM  87 
JMGM  88 
JMCM  89 
JMGM  90 
JMGM  91 
JMGM  92 
JMGM  93 
JMGM  94 
JMGM  55 
JMGM  95 
JMCM  97 
JMGM  98 
JMGM  99 
JMGM  100 
JMGM  101 
JMGM  102 
JMGM  103 
JMGM  104 
JMGM  105 
JMGM  10B 
jrSM  107 
JMGM  103 
JMGM  109 
JMGM  110 
JMGM  111 

JMGM  112 

, JMGM  113 

JMGM  114 
GIUENJMGM  115 
JMGM  115 
JMGM  117 
JMGM  118 
JMGM  119 
JMGM  120 
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c 
c. 
c. 
c 

c 
c 


c 

c 
c 
c 
c 
c 
c 
c 
c 
c 

c 
c 
c 

c 

c 


c 

c 


c 

c 


c 

c 
c 
c 
c 
c 


c 
c 
c 


c 

c 

c 


CITER=2.AITER=0.IZ=1    FIND  BETA(I)  S  FROM  REGRESSION  EQUA. 

300  RTP=UPARMSQ)  +  (UPARMS(2)-1)«3 
ATT=0.0 
TUAR=0.0 
SX=0.0 
SXX=0.0 

DO  18  J=1.0MAX 
IF(Z(2.J'.EO.0.0)GO  TO  18 

COMPUTE  RAUE,   AUERAGE  FROM  REGRESSION  EQUATION 

AUE  =  ZONAL  AUERAGE  TRAUEL  TIME. 

AUEO  =  ZONAL  AUERAGE  TRAUEL  TIME  AT  BETA  =0.0 

UAR   =  ZONAL  TRAUEL  TIME  UARIANCE  AT  CALIBRATED  BETA  UALUE. 

UARO  =  ZONAL  TRAUEL  TIME  UARIANCE  AT  BETA  =0.0 

TERM  =  AUERAGE  ORIGIN  ZONE  TERMINAL  TIME. 

GO  TO(  1.2, 3, 4, 5. 6, 7. 8, 9. 10, 11. 125.RTP 

COMBINED  REGRESSION  EQUATIONS-LAFAYETTE,  ANDERSON.  AND  MUNCIE  IND. 

HBU   R»»2  =  0.8S5 

1  RAUE=-0. 4327+1. 2761*Z(4.J)+0.S921«(Z(10.J)-Z(4,J» 

GO  TO  19 

HBO   R»*2  =  0.877 

2  RAUE=-0.7499+1.3548*Z(4.J)+0.7557»(Z(10.JW(4.J» 

GO  TO  19 

NHB   R**2  =  0 . 735 

3  RAUE=0. 5245+1. 3416«Z( 4.  J)+0.G073»(ZQ0.J)-Z(4,jn 

GO  TO  19 

LAFAYETTE,  INDIANA. 

HBU   R»*2  =  .893    AREAUIDE  UAR=22.99:   UAR0=24.99 
AUE=10.G0:   AUE0=12.32:   UAR=12.17:   UARC=14.00;   TERM=1.08 

4  RAUE=-1.0SS1+1.0934*Z(4,J)+0.3328»(Z(10,J)-Z(4,J)) 
GO  TO  19 

HBO   R»*2  =  .893    AREAUIDE  UAR=13.B3:   UAR0=25.5^ 
AUE=8.55:   AUE0=11.78:   UAR=10.78S   UAR0=15.63i   TERM=1.11 

5  RAUE=-1. 6791+1. 145S»Z(4,J)+0.8397»CZC10.J)-Z(4,J» 
GO  TO  19 

NHB   R**2  =  .799    AREAUIDE  UAR=14.09:   UAR0=17.75 
AUE=9.12;   AUE0=11.42;   UAR=9.94:   UAR0=12.9l:   TERM=1.93 
S  RAUE=0.2131+1.2823*Z(4,J)+0.6744»(Z(10.J)-Z(4,J)) 
GO  TO  19 


C 
C 
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ANDERSON,  INDIANA 
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JMGM 
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13S 
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JMGM 
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JMGM 
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JMGM 
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JMG^ 
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JMGM 

1^2 

JMGM 

1-3 

J.-,3M 
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1^5 
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,JMGM 

147 

JMGM 
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jMGM 
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JMGM 

15-1 
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158 
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JMGM 
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JMGM 
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JMGM 
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JMGM 
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JMGM 
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JMGM 

1SS 
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JMGM 
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jMGM 
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JMGM 
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173 

JMGM 
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JMGM 
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JMGM 

177 

JMGM 

:" 

JMGM 
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c 
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c 
c 
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c 

c 
c 

c 
c 
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c 

c 
c 


c 
c 
c 


c 
c 
c 


c 

c. 

c. 

c 
c 
c 


HBW   R»«2  =  .876    AREAUIDE  UAR=9.61;   UAR0=9.51 
AUE=11.43;   AUE0=12.10:   UAR=6.53:   U0=7.12:   TERM=2.38 

7  RAUE=-0.4650+1.1306»Z(4.J)+0.9467»(Z(10.J)-Z(4,J)) 
GO  TO  19 

HBO   R**2  =  .664    AREAUIDE  UAR=9.51:   UAR0=11.10 
AUE=10.04:   AUE0=12.30:   UAR=7.84:   UAR0=5.20:   TERM=2.37 

8  RAUE=1.0G11+1.0640»Z(4,J)+0.6504»(Z(10, J)-Z(4,J» 
GO  TO  19 

NHB   R»»2=  .210    AREAUIDE  UAR=7.82;   UAR0=9.20 
AUE=9.50:   AUE0=11.72;   UAR=G.35:   L'AR0=7.12;   TERM=2.G4 

9  RAUE=274556+1.0868»Z(4.J)+0.4280«(Z(10.J)-Z(4.J)) 
GO  TO  19 


MUNCIE.  INDIANA 

HEW   R»»2  =  .772    AREAUIDE  UAR=11.16;   UAR0=11.51 
AUE=9.07;   AUE0=9.82!   UAR=7.06;   UAR0=7.11:   TERM=1.33 

10  RAUE=0.8G18+0.7448»Z(4.J^+0.8491«(2(10,J)-2(4,J)) 
GO  TO  19 

HBO   R**2  =  .749    AREAUIDE  UAR=9.27;   UAR0=11.83 
AUE=7.48:   AUE0=9.58:   UAR=6.40;   UAR0=7.79:   TERM=1.34 

11  RAUE=0. 1603+1. 0595«Z(4,J)+0.7162»(Z( 10, J)-ZC4,J» 
GO  TO  19 

NHB   R*«2  =  .506    AREAUIDE  UAR=7.89:   UAR0=9.28 
AUE=6.93;   AUE0=8.95:   UAR=5.65;   UAR0=G.G5:   TERM=1.49 

12  RAUE=0. 7980+1. 4192«Z(4. J)+0.5384«(ZC 10.  J)-Z(4. J)) 

19  CONTINUE 

COMPUTE  BETA,  LINEAR  INTERPOLATION. 

Z(G,J)=(PBETA(TP)/'(Z(12,J)-Z(10,J)))»(RAUE-Z(10.J)) 
IF(Z(12,J).GT.Z(10.J))Z(G,J)=»0.0 
IFCZtG. J).GT.0.5)Z(6.J)=0.S 
18  CONTINUE 
GO  TO  99 


AITER=0,IZ=1 

400  DO  401  J=1.0MAX 
Z(8.J)=Z(3.J) 

401  2(9,J)=0.0 
AITER=1 


INITIALIZE  BEFORE  ATTRACTION  ITERATION  BEGINS. 


C 
C 

c 
c 
c 


BEGIN  ATTRACTION  ITERATION. 

SUMJ  A(J)»F(I.J) 

500  Z(7,IZ)=0.0 

DO  501  J=1.0MAX 
IF(IZ.E0.J)INT(J.1)=Z(5.J) 
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223 

JMGM 

22- 

JMGM 

225 

J^GM 

225 

JMGM 

227 

J"GM 

223 

JMGM 

229 

JMGM 

230 

JMGM 

231 

JMGM 

232 

JMGM 

233 

JMGM 

234 

J'-GM 

235 

JMGM 

235 

JMGM 

237 

JMGM 

233 

JMGM 

239 

JMGM 

240 

40 


501 


C 

C 

c 


502 


C 

c 
c 
c 


c 
c 
c 


700 
701 


C 
C 
C 

C 

c 
c 


703 


c 
c 
c 

c 
c 
c 


c. 

c 

c 

c 
c 
c 


702 


IF(IW(J. l).EQ.0.0)GO  TO  501 

Z(14p J)  =  1.0 

ADD  TERMINAL  TIMES  TO  SKIM  TREE  TRAUEL  TIME. 

TiME=INT( J. 1)+Z(4, I2)+Z(4,  J) 

IF(TIi1E.GT.99>TIME=93 

IF(CITER.E0.1)Z(14,J)=FCTF(TIME) 

IF(CITER.EQ.a)Z(14tJ)=EXPC-Z(6.IZ)«TIME) 

SAF  =  SAF  +  A  •  F 

Z(7.IZ)=Z(7,IZ)+Z(8,J)*Z(14,J) 

CONTINUE 

SUM  COMPUTED  ATTRACTIONS 

DO  502  J=1.0MA!< 
IF(INT(J.l).EQ.0.0)GO  TO  502 
SA=SA+P»A»F/  SAF 
Z(3.J)=Z(3.J)+Z(2.IZ)*Z(8,J)»Z(14,J)/'Z(7,IZ) 

CONTINUE 

IF  IZ.EO.OMAX.   TEST  FOR  CLOSURE 
IF  IZ.LT.OMAX.   CONTINUE  ITERATION 

IF(IZ.EQ.OMAX)GO  TO  700 
RETURN  1 

COMPUTE  ERROR  FOR  CLOSURE  TEST 

SUM=0.0 

DO  701  J=l,OMAX 

SUM=SUM+ABS(Z(9.J)-Z(3.J)) 

SUM=SUM/TOT 

ONE  PERCENT  CLOSURE  TEST 

IFCSUM.LT.0.0DGO  TO  702 

CLOSURE  TEST  FAILED,   COMPUTE  NEW  ATTRACTIONS. 

DC  703  J=1.0MAX 
IF(Z(8.J).EQ.0.0)GO  TO  703 
Z(8,J)=Z(3,J)»Z(8,J)/Z(9.J) 
Z(9.J)=0.0 
RETURN  1 

CLOSURE  TEST  PASSED.  SET  AITER=2  FOR  AUE.  AND  UAR.  ITERATION. 

AITER=2 

IF  CITER=2,  SET  ITNG=ITER-1  FOR  LAST  AUE.  UAR.  ITERATION. 

IF(CITER.EQ.2)ITN0=ITER-1 
RETURN  1 


AUERAGE  AND  UARIANCE  COMPUTATION  . 
INITIALIZE  ZONAL  AUE.  AND  UAR. 


JMGM  2*1 

JMGM  2-'? 

JMCM  343 

2^4 

JMGM  245 
J.1GM  24  S 
247 
JMGM  248 
uMGM  249 
JMGM  250 
JMGM  251 
JMGM  252 
JMGM  253 
JMGM  254 
JMGM  255 
JMGM  256 
JMGM  257 
JMGM  253 
JMGM  259 
JMGM  260 
JMGM  261 
JMGM  252 
JXM  253 
JMGM  264 
JMGM  255 
JMGM  266 
JMGM  267 
JMGM  266 
JMGM  259 
JMCM  270 
JMGM  271 
JMGM  272 
JMGM  273 
JMGM  274 
JMGM  275 
JMGM  275 
JMGM  277 
JMGM  278 
JMGM  279 
JMGM  280 
JMGM  281 
JMGM  2S2 
JMGM  283 
.'3M  234 
JfIGM  235 
JMGM  2S6 
JMGM  287 
JMGM  28E 
JMGM  293 
JMGM  230 
JrSGM  291 
JMGM  292 
_~:~  2=3 
..JMGM  23d 
..JMGM  235 
JMGM  2S5 
JMGM  297 
JMGM  29S 
JMGM  299 
JMGM  300 
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C 
C 

c 
c 
c 


GOO  AUE=O.0 
UAR=0.0 
DO  539  J=l. ZONES 

5"  lFC((cifERTEQ.O).OR.(CITER.EO.a)).AND.(IZ.E0.1))CO  TO  598 

GO  TO  597 
593  ATT=0.0 
TUAR=0.0 
SX=0.0 
SXX-0.0 

BEGIN  AUE.  AND  UAR.  COMPUTATION 

IF  ORIGIN  ZONE  PRODUCTIONS.  SKIP  AUE.  AND  UAR. 


C 
C 

C 


c 

c 
c 

c 
c 
c 


c 
c 
c 
c 

c 
c 
c 

c 

c 
c 

c 
c 
c 

c 
c 
c 


597  IF(Z(2.IZ).LE.1.0)GO  TO  605 
DO  SCI  J=1.0MAX 

CONFUTE  CELL  INTERCHANGE.  I.E.  TRIPS. 

IF(IZ.EQ.J)INT(J.1)=ZC5»J) 
IF(INT(J.l).EQ.0.0)GO  TO  601 

Z(14,J)=1 
TIME=INT(J.1)+Z(4, IZ)+Z(4.J) 

IF(TIME.CT.99)TIME=99 

IF(CITER.E0.1)Z(14.J)=FCTR(TIME) 

IFCCITER.EQ.2)ZC14,J)=EXP(-Z(G.IZ)»TIME) 

TRIP  =  P  *  A  •  F  /  SAF 
TRIP=Z(2.IZ)»Z(8.J)*Z(14.J)/-Z(7,IZ) 

SAUE  TRIP  IN  TROUT  IF  LAST  ITERATION,  I.E.  CITER=2 

IF( (CITER.EQ.2) .AND. (AITER.EQ.2) )TROUT( J, 1J-TRIP 

SUM  ZONAL  TRAUEL  TIMES. 

AUE=AUE+TRIP«TIME 
BOl  CONTINUE 

ADD  ZONAL  CONTRIBUTION  TO  AREAWIDE  AUERAGE  TRAUEL  TIME 
IFC(CITER.EQ.O). OR. (CITER.EQ.2)) 

IF((CITER.EQ.0).OR.(CITER.EQ.2))ATT=ATT+AUE/TOT 

COMPUTE  ZONAL  AUERAGE  TRAUEL  TIME 

AUE=AUE/Z(2.IZ) 

SAUE  AUERAGE  TRAUEL  TIME  AT  BETA=0.0  IF  CITER=0 

IF(CITER.EO.0)Z(10,IZ)=AUE 

SAUE  AUERAGE  TRAUEL  TIME  AT  BETA=PBETA(TP)  OR  FINAL  UALUE. 

IF(CITER.NE.0)Z(12»IZ)=AUE 

COMPUTE  ZONAL  UARIANCE  AND  SAUE  DATA  FOR  AREAWIDE  UARIANCE 

DO  602  J=1.0MAX 


JMGM  301 
JMGM  302 
JMGM 

JMGM  304 
JMGM  305 
JMGM  30E 
JMGM  307 
JiIGM  303 
JMGM  309 
JMGM  310 
JMGM  311 
JMGM  212 
JMGM  313 
JMGM  314 
JMGM  315 
JMGM  315 
JMGM  317 
JMGM  318 
JMGM  319 
JMGM  320 
JMGM  321 
JMGM  322 
JMGM  323 
JMGM  324 
JMGM  325 
JMGM  326 
JMGM  327 
JMGM  32S 
JMGM  323 
JMGM  330 
JMGM  331 
JMGM  332 
JMGM  333 
JMGM  334 
JMGM  335 
JMGM  33G 
JMGM  337 
JMGM  338 
JMGM  339 
JMGM  340 
JMGM  341 
JMGM  342 
JMGM  343 
JMGM  344 
JMGM  345 
JMGM  34G 
JMGM  347 
JMGM  348 
JMGM  343 
JMGM  350 
JMGM  351 
JMGM  352 
JMGM  353 
JMGM  354 
JMGM  355 
JMGM  356 
JriGM  357 
JMGM  353 
JMGM  35S 
JMGM  3G0 
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COMPUTE  CELL   INTERCHANGE.    I.E.   TRIP 

IF(INT(J, l).EQ.0.0)GO  TO  SOS 
TRIP=2(2. IZ)«Z(8.  J)«Z(14,  JVZC7, IZ) 

ADD  SUN  OF  SQUARES  FOR  ZONAL  UARIANCE 

UAR=UAR+TRIP»(TIME-AUE)«»2 

ADD  SUN  OF  TINE  AND  SUN  OF  TINE  SQUARED  FOR  AREAWIDE  UARIANCE 
IF(CCITER.EQ.0).0R.(CITER.EQ.2)) 

IF(INT(J.l).EO.0)GO  TO  G02 
TINE=INT(  J.  1 )+Z(4,  IZKZC4,  J) 
IF(TIME.GT.99)TIME=99 

IF((CITER.EQ.0).0R.(CITER.EQ.2))SX=SX+TRIP»TIME 
IF((CITER.EQ.0).0R.(CITER.EQ.2'»)SXX=SXX+TRIP»TIME»«2 
S02  CONTINUE 

SAUE  ZCNAL  UARIANCE  AT  BETA=0.0  IF  CITER=0. 

IF(CITER.EQ.0)Z(11.IZ)=UAR/(Z(2. IZ)-1.0) 

SAUE  ZONAL  UARIANCE  AT  BETA=PBETA(TP)  OR  FINAL  UALUE. 

IF(CITER.NE.0)ZC13.IZ)=UAR/(Z(2.IZ)-1.0) 


C 

c 
c 

c 

C 

C 


IF  IZ=ONAX  INCREMENT  CITER  BY  1  AND  SET  AITER=0.   STMT  S03 

IF(IZ.EQ.OMAX)GO  TO  G03 

IF  IZ.LT.OMAX.   CONTINUE  ITERATION. 

RETURN  1 


G03  CITER=CITER+1 

IF(CITER.EQ.1)URITE(6.SCG)ATT.TUAR 

EOG  FORNATC/10X.10HBETA   =   0.  (V10X, 22HAUERAGE  TRAUEL   TIME  =   .F10.3/ 

H0X.23HTRAUEL   TIME  UARIANCE  =   .F9.3/0 

IF(CITER.EG.3)URITECS.604)ATT,TUAR 

G04  FORMAT(/10X.22HAUERAGE  TRAUEL  TIME  =  , F10.3/10X.23HTRAUEL  TIME  UAR-TGM 

1IANCE  =  .F9.3/)  -MGM 

AITER=0  JMGM 

RETURN  1  -"-•' 
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